A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

11.02.2: Mkn 501;12.04.2;13.09.2;13.07.2;13.25.2 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



O ■ 

O ' 

o : 

(N ■ 

<\ 

<n : 



> ■ 

in : 
in . 
m ■ 
■ 

o : 
o . 
o ■ 

^ ' 

Or 

6 ■ 
a: 

c/j . 



Constraints on the Cosmic Infra-Red Background based on 
BeppoSAX and CAT spectra of Mkn 501 

J. Guy 1 , C. Renault 1 , F.A. Aharonian 2 , M. Rivoal 1 , J.-R Tavernet 1 

1 LPNHE, CNRS-IN2P3 Universites Paris VI- VII, 4 place Jussieu, F-75252 Paris Cedex 05, France 

2 Max Planck Institut fur Kernphysik, Postfach 103980, D-69029 Heidelberg, Germany 

Received; Accepted 



Abstract. The TeV and X-ray data obtained by the 
imaging Cherenkov telescope CAT and X-ray satellite 
BeppoSAX during the remarkable flare of Mkn 501 in 
April 16, 1997 are used to constrain the flux of the Cos- 
mic Infrared Background (CIB) using different CIB mod- 
els. We show that a non-negligible absorption of 7-rays 
due to the CIB could take place already in the low-energy 
(sub- TeV) domain of the spectrum of Mkn 501. This im- 
plies that the data of the low-energy threshold CAT tele- 
scope contain very important information about the CIB 
at short wavelengths, 0.4 fim< A < 3. fim. The analy- 
sis of almost simultaneous spectroscopic measurements 
of Mkn 501 in a high state by CAT and BeppoSAX in 
the framework of the standard homogeneous Synchrotron- 
Sclf-Compton (SSC) framework model leads to the con- 
clusion that the density of the near-infrared background 
with typical "starlight spectrum" around 1 /im should 
be between 5 and 35nWm _2 sr _1 (99% CL), with most 
likely value around 20 nWm _2 sr _1 . Also we argue that 
the CAT 7-ray data alone allow rather robust upper lim- 
its on the CIB, \F\ < 60 nWm^sr 1 at 1 /im, tak- 
ing into account that for any reasonable scenario of 7- 
ray production the differential intrinsic spectrum of 7- 
rays hardly could be flatter than dN/dE oc E^ 1 . This 
estimate agrees, within statistical and systematic uncer- 
tainties, with recent reports about tentative detections of 
the CIB at 2.2 and 3.5 /im by the Diffuse Infrared Back- 
ground Experiment (DIRBE), as well as with the mea- 
surements of the background radiation at optical wave- 
lengths from absolute photometry. The high flux of CIB 
at < few /im wavelengths implies a significant distortion 
of the shape of the initial (source) spectrum of 7-rays from 
Mkn 501 at sub- TeV energies. The "reconstructed" intrin- 
sic 7-ray spectrum shows a distinct peak in the Spectral 
Energy Distribution (SED) around 2 TeV with a flux by 
a factor of 3 higher than the measured flux. The energy 
spectrum of gamma radiation from both sides of the peak 
has power-law behavior with photon index a ~ 1.5 be- 



low 2 TeV, and a ~ 2.5 above 2 TeV. This agrees with 
predictions of SSC model. We also discuss the impact of 
the intergalactic absorption effect in derivation of the SSC 
parameters for the jet in Mkn 501. 
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1. Introduction 

The Cosmic Infrared Background (CIB) is contributed 
mostly by the red-shifted 'stellar' and 'dust' radiation 
components, and therefore carries vital cosmological in- 
formation about the epoch of galaxy formation and their 
evolution in time. The derivation of information about 
CIB from direct measurements is a hard task which re- 
quires an effective removal of heavy contamination caused 
by foregrounds of different origin, in particular by the zo- 
diacal light, stellar and interstellar emission of our Galaxy, 
etc (Arendt et al. 1998; Kelsall et al. 1998). The conclu- 
sions of this approach in the near-infrared domain are to 
a large extent model-dependent because they are gener- 
ally based on comprehensive modeling of the foregrounds; 
the far-infrared background is determined with a better 
accuracy as it dominates the foreground emission (Hauser 
et al. 1998). 

Very High Energy (VHE) 7-ray astronomy provides an 
independent and complementary approach for the study 
of the CIB. The idea is simple, and based on the detec- 
tion of absorption features in the 7-ray spectra of distant 
extragalactic objects caused by interactions of VHE 7- 
rays with the CIB photons in their way from a source to 
the observer (Nikishov 1962; Gould and Schreder 1967, 
Stecker et al. 1992). The recent detections of 7-rays from 
two BL Lac objects, Mkn 421 and Mkn 501, with spectra 
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extending up to 10 TeV and beyond, open an interest- 
ing path for the realization of this exciting cosmological 
aspect of VHE gamma-ray astronomy. Obviously, the suc- 
cess of the '7-astronomical' approach essentially depends 
on two crucial conditions: (i) accurate 7-ray spectrometry, 
and (ii) good understanding of the intrinsic (source) spec- 
tra of TeV 7-rays. 

The observations of Mkn 501, the second closest X-ray 
selected BL Lac object with a redshift z ~ 0.034, during 
its extraordinary outburst in 1997 yielded unique data, 
which initiated important theoretical studies of the phys- 
ical conditions in the relativistic jets of BL Lac objects 
(see e.g. Tavecchio et al. 1998 (hereafter TMG); Kirk & 
Mastichiadis 1999, Hillas 1999; Coppi & Aharonian 1999a, 
Bcdnarek & Protheroe 1999, Krawczynski et al. 1999), as 
well as interesting efforts to set meaningful upper limits 
on the CIB flux (see Biller et al. 1998, Stanev & Frances- 
chini 1998, Barrau 1998, Stcckcr & De Jager 1998, Stecker 
1999, Aharonian et al. 1999, Coppi & Aharonian 1999b; 
Konopelko et al. 1999). 

During the 1997 outburst, lasted several months, 
Mkn 501 showed dramatic variations of fluxes both in 
X-rays (BeppoSAX: Pian et al. 1998, RXTE: Lamer & 
Wagner 1998) and TeV 7-rays (Whipple: Catanese et al. 
1997, Samuelson et al. 1998; HEGRA: Aharonian et al. 
1997, 1999; Telescope Array: Hayashida et al. 1998; CAT: 
Djannati-Ata'i et al. 1999). More importantly, the high 
TeV fluxes allowed monitoring of the energy spectrum of 
Mkn 501 on a diurnal basis, especially during very strong 
flares. On several occasions truly simultaneous observa- 
tions of Mkn 501 were available in TeV and X-ray band 
(Krawczynski et al. 1999). A special interest presents 
the famous April 16, 1997 flare which was observed by 
BeppoSAX (Pian et al. 1998) and low-energy threshold 
(<~ 300 GeV) Whipple and CAT atmospheric Cherenkov 
telescopes (Catanese et al. 1997; Djannati-Ata'i et al. 
1999). Because of very large fluxes, the energy spectrum 
of the flare was obtained with good accuracy in broad dy- 
namical ranges in both X-ray (0.1-100 keV) and TeV 7- 
ray (0.3-10 TeV) regimes. Within today's most popular, 
the so-called synchrotron-self-Compton (SSC) model of 
non-thermal high energy radiation of Mkn 501, the data of 
April 16, 1997 flare allow to compute the intrinsic source 
spectrum of 7-rays. This information, coupled with al- 
most simultaneous spectral measurements of 7-radiation 
by CAT, provides a good opportunity to analyze the in- 
tergalactic absorption signature in the observed VHE 7- 
ray spectrum. 

Since the optical depth r 77 of 7-rays in the intergalac- 
tic medium increases with energy (for any reasonable 
spectral shape of the CIB), most stringent constraints 
on the CIB come from the very energetic tail of the 7- 
ray spectrum of Mkn 501. For E > 10 TeV the 7 — 7 ab- 
sorption is dominated by the 'dust' component of CIB at 
far infrared (FIR) wavelengths, A > 10 /mi. Therefore, it 
is generally believed that a deep probe of CIB at optical 



and near-infrared (NIR) wavelengths is contingent only 
upon the discovery of more distant (e.g. with z > 0.1) 
VHE 7-ray sources. In this paper, however, we show that 
the accurate measurements of the spectrum of Mkn 501 
by CAT at sub- TeV energies already provide meaningful 
upper limits on the CIB flux at wavelengths between 0.4 
and a few ^m. Moreover, the analysis of the X-ray and 
7-ray spectra of the April 16, 1997 flare within the ho- 
mogeneous SSC model allows rather conclusive estimates 
of the CIB flux at such short wavelengths. The CIB flux 
\F\ ^5 — 35 nW m _2 sr _1 at A <~ 1 /mi gives a reasonable 
slope, vF v oc E - 5 , in the 'reconstructed' Spectral Energy 
Distribution (SED) of 7-rays at low energies - as expected 
within the framework of the SSC model. The significant 
intcrgalactic absorption of sub- TeV 7-rays leads to both 
the shift of the position (E ~ 2 TeV ) and increase of the 
flux [yF v ~ 10~ 9 erg cm~ 2 s _1 ) of the so-called Compton 
peak in the 7-ray spectrum. An analytical approach has 
been recently proposed by (TMG) for derivation of con- 
straints on the jet parameters of TeV blazars. With these 
revised spectral parameters, new results in a well defined 
self-consistent SSC parameter-space for the jet of Mkn 501 
in the high state are obtained. 

2. Intergalactic absorption of VHE gamma rays 

If we ignore the appearance of second-generation 7- 
ray photons in the source direction 1 , the extinction of 
gamma radiation is reduced to a simple absorption effect, 
described by an energy dependent optical depth: 

J obs {E) = ME) cxp [-t 7J (E)\ 

where Jo and J b s are the intrinsic (source) and detected 
7-ray fluxes, respectively. In order to calculate the optical 
depth for 7-rays from several hundred GeV to 20 TeV (the 
typical energy range of current Cherenkov telescopes), one 
needs to know the distance to the source d = cz/Hq, 
where H n is the Hubble constant, and the flux of the 
diffuse background in a rather broad interval from sub- 
micron to 100 /mi wavelengths (for close sources, z <C 1, 
the CIB flux is independent of z). Note that if, for a given 
redshift z, the uncertainty in the estimate of the distance 
to the source does not exceed a factor of 1.5 (the current 
estimate of the Hubble constant Hq is between 50 and 

1 This is a good approximation unless the intergalactic mag- 
netic field B < 1(T 16 G. In the case of such extremely small 
intergalactic magnetic fields, the new generation (cascade) 7- 
rays, produced by secondary electrons and positrons via inverse 
Compton scattering on 2.7 K microwave background, remain 
in the field of view of the detector and arrive almost simulta- 
neously with the primary 7-rays. If so, the cascade spectrum 
could dominate over the primary 7-ray spectrum, and therefore 
instead of the absorption features in the primary spectrum, we 
should detect a standard power-law (oc E~ ) cascade spec- 
trum with exponential cutoff at E* , determined from the con- 
dition t{E*) = 1. 
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75 km/s/Mpc), the uncertainty connected with the CIB 
flux (directly measured and/or predicted by different cos- 
mological models) is significantly larger. However, since 
in the method described below we intend to obtain inde- 
pendently the density of the CIB using the 7-ray data, we 
need to know the spectral shape (or possible shapes) of 
the CIB, rather than the absolute value of the CIB, in a 
broad band of wavelengths from 0.3 /im to 100 /im. 

2.1. CIB models 

The calculations of CIB present a serious theoretical chal- 
lenge because they require a complex treatment of number 
of cosmological assumptions and key physical processes 
of galaxy formation (see e.g. MacMinn & Primack 1996; 
Salamon & Stecker 1998, Madau 1999, Primack et al. 
1999). It is important for our further discussion that most 
of cosmological models give rather similar shapes of the 
basic, the 'stellar' and 'dust' components of radiation - 
two distinct bumps at 1-2 /im and 100-200 /j,m (see e.g. 
Dwekctal. 1998; Primack et al. 1999). The shape (the 
precise position, depth and width) of the mid-infrared 'val- 
ley' contains more uncertainties because of the lack of ad- 
equate information necessary for modeling the radiation 
associated with warm dust component. 

In this paper, we use a sample of recent models pro- 
viding information from roughly 0.1 to 300 /im. The pre- 
dictions of two such models, LCDM (Lambda Cold Dark 
Matter) and CHDM (Cold and Hot Dark Matter) sug- 
gested by Primack et al. (1999) for the density of CIB 
is shown in Fig.l (top panel). For comparison we show 
also the maximum and minimum CIB fluxes derived phe- 
nomcnologically by Malkan & Stecker (1998) for wave- 
lengths A > 3.5/im. However, in this paper we are pri- 
marily interested in absorption of low energy, E < 2 TeV 
7-rays, which effectively interact with NIR background 
photons. Therefore below we will not use the CIB models 
of Malkan & Stecker (1998). 

In the bottom panel we present another three CIB 
models suggested by Dwek et al. (1998) which differ from 
each other by cosmic star formation histories. The PFI and 
PFC models are based on different chemical evolution sce- 
narios directly related to the star formation rate based on 
the Pei & Fall (1995) calculations. The ED model (Extra- 
galactic DIRBE results) is derived from UV and optical 
observations and reproduce the COBE fluxes (for details 
see Dwek et al. 1998). 

The flux dispersion between the curves in Fig.l shows 
that neither the shape nor the intensity can be precisely 
predicted by current CIB models. Actually in our ap- 
proach we intend to "measure" the CIB density ourselves. 
Therefore we are interested in the shape rather than in 
the absolute flux predicted by CIB models. Fortunately, 
in the most informative parts of the gamma-ray spectrum, 
namely below 1 TeV and above 10 TeV, the absorption of 
7-rays from relatively nearby BL Lac object Mkn 501 is 



contributed mainly by two well separated parts of the CIB 
spectrum - A < few /jm and A > 10 /im, respectively. This 
allows us to introduce a scaling factor SF, which does 
not change the shape, but varies the absolute flux of CIB 
within reasonable limits, and thus to derive the "best-fit" 
values or upper limits on SF for each specific CIB model. 

2.2. 7 — 7 pair production 

A 7-ray photon with energy E penetrating through an 
isotropic field of photons, can interact with any back- 
ground photon of energy e > t t h — {m e c 2 ) 2 / E ~ 
0.26 (E/l TeV) -1 eV. Since the cross-section of the pair 
production 77 — ► e + e~ peaks at e max — ^eth with 
cr 77 ~ 10~ 25 cm 2 (sec e.g. Herterich 1974), for a large class 
of relatively flat and smooth (e.g. power-law) spectra of 
the field photons, approximately half of the optical depth 
r 77 is contributed by a rather narrow band of the back- 
ground radiation within e max ± l/2e max . For a more real- 
istic shape of the CIB photons with two distinct bumps at 
NIR and FIR, the relative contributions of different parts 
of the spectrum of CIB do not follow this simple relation. 
In Fig. 1, we show the energy intervals from the thresh- 
old tth of CIB photons contributing to the 50% and 90% 
of the optical depth of intcrgalactic absorption calculated 
for the LCDM and the ED type CIB spectra for 3 ener- 
gies of primary 7-rays: E =600 GeV, 4 TeV, and 17 TeV. 
It is seen that 50% of absorption of 600 GeV 7-rays is 
caused by a narrow band of CIB between 1-3 /im, while 
the remaining contribution comes essentially from the 0.4- 
1 /itm band for the LCDM model. The thresholds depends 
only on the TeV photon energy, so they are independent 
of the CIB model. The bands are significantly narrower 
for the ED model because both NIR and FIR parts are 
more "peaked". 

The absorption of multi-TeV 7-rays is equally con- 
tributed both from mid- and far-IR parts of the spectrum. 
So conclusions based on the absorption of these energetic 
7-rays depend strongly on the model of the CIB, first of 
all on the depth of the mid-IR 'valley' of the CIB. It is im- 
portant to note that conclusions concerning the change of 
the spectral shape of sub- TeV 7-rays due to the CIB ab- 
sorption is independent of the density of mid- and far-IR 
photons, these photons being beyond the 7 — 7 interaction 
threshold. 

The optical depth is a function of the 7 photon energy 
E defined as, for a source at a redshift z s « 1 : 

r( £ ) = |£/ L^pl d{cos0) n(e)a(E,e,9)de 

J — 1 - J e £ 

where 9 is the angle between both photons, e the IR 
photon energy and a(E, e, 9) is the pair-production cross- 
section. 

In Fig. 2 we present the portion (percentage) of ab- 
sorbed 7-rays, P = 1 — exp[— t 77 (-E)], from Mkn 501 for 
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Fig. 1. Predictions for the CIB fluxes by different cosmologi- 
cal models, a: the LCDM and CHDM models of Primack et al. 
(1999) are shown by thick solid and dashed lines, respectively. 
The maximum and minimum CIB fluxes from the model of 
Malkan and Stecker (1998) are shown by thin solid and dashed 
lines, respectively. The arrows show the wavelength intervals 
(from the threshold) of CIB photons with 50% and 90% contri- 
butions to the absorption of TeV 7-rays with energies 600 GeV, 
4 TeV, and 17 TeV computed for the LCDM distribution, b: 
three models discussed by Dwek et al. (1998) are drawn with 
solid line (PFI model) , dashed (PFC model) and grey line (ED 
model). The arrows show the wavelength intervals computed 
for the ED distribution. The thresholds indicated by the left 
side of the arrows is independent of the CIB model. 



the LCDM and ED models assuming two different scal- 
ing factors: SF — 1 and 2.5. Hereafter all calculations are 
performed for the Hubble constant Hq — 60km/sMpc. 

3. Constraints on CIB 

Given the lack of reliable information on the intrinsic 7- 
ray spectrum, an upper limit on the CIB flux could be 
derived by formulation of a priori, but astrophysically 
meaningful requirement on the shape of the spectrum of 



Fig. 2. Absorption in percents of 7-rays from Mkn 501 in func- 
tion of the gamma energy calculated for two different scaling 
factors (1 and 2.5) and two models (LCDM and ED). 



7-rays produced in the source. In particular, we may re- 
quire that within any reasonable model of 7-ray produc- 
tion, and for a given observed 7-ray spectrum J b s {E), the 
source spectrum, Jo(E) — J b s (E) exp[r 77 (E)], should not 
contain, at any energy E, a strongly (e.g. exponentially) 
rising feature. 

In Fig. 3 we show the spectrum of Mkn 501 as 
measured by the CAT telescope during the strong 
April 16, 1997 flare in the energy region from 300 GeV to 
12 TeV (statistical errors only). Unfortunately, because of 
bad weather this remarkable flare could not be observed 
by the HEGRA telescope system. However, the HEGRA 
observations of Mkn 501 revealed that, despite dramatic 
flux variations in time-scales < 1 day, the shape of the 
energy spectrum above 1 TeV remained essentially stable 
throughout the entire state of source high activity in 1997. 
So we show the HEGRA 'time-averaged' spectrum (Aha- 
ronian et al. 1999) normalized to the CAT April 16, 1997 
spectrum at E = 1 TeV with the re-scaling factor of «2.2. 
At TeV energies the agreement between the CAT and 
HEGRA spectra is quite impressive. Below 1 TeV we show 
only the CAT spectral points since in this energy region 
both the statistical and systematic errors of the data ob- 
tained close to the energy threshold of the HEGRA tele- 
scope system are rather large. Moreover, the inclusion of 
the HEGRA 'time-averaged' spectral points at sub- TeV 
energies in Fig. 3 cannot be justified because of reported 
noticeable variation of the spectrum of Mkn 501 at such 
low energies in the time-scales < 1 day, especially during 
strong flares (Djannati-Ata'i et al. 1999). 

3.1. Limit on the CIB at FIR 

More than half of the intergalactic absorption of high- 
est energy photon above 10 TeV detected from Mkn 501 
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Fig. 3. Spectral energy distribution of 7-rays from Mkn 501 
as detected by CAT during the April 16, 1997 flare (filled 
squares). The spectral points of the HEGRA 'time-averaged' 
spectrum normalized to the CAT flux at 1 TeV are also pre- 
sented (open squares). Statistical errors only are shown. Three 
intrinsic 7-ray spectra are computed using the LCDM model 
and scaling factors of 1, 2.5 and 4. 



tradict this 'reasonable-source-spectrum' criterion. Actu- 
ally, the acceptable range of SF at 68% C.L. extends to 
SF = 3.4 (see Tab. 1). Note that as the above formulated 
i/i^(17TeV)/i'F 1/ (4TeV) < 1 criterion could be treated as 
a tight condition, the results presented in Tab.l are based 
on a rather relaxed, 68% C.L. statistical requirement. 

The upper limits on SF are quite sensitive to the 
exact shape of the CIB spectrum because the absorp- 
tion of 17 TeV 7-rays is contributed mostly by the back- 
ground photons from the rising branch of the CIB spec- 
trum (A > 10 /im, see fig. 1), therefore the optical depth 
increases strongly with energy of 7-rays. 

The results are summarized in Tab. 1. The derived val- 
ues of SF ma x give upper limits on the absolute flux of CIB 
from 20 to 80 ^m shown in Fig. 5. Note that the CHDM, 
LCDM, and ED models have rather similar slopes between 
20 and 80 fim (see Fig. 1). Therefore the upper limits on 
CIB corresponding to these three models are very close. 
The same is true for the PFI and PFC models. For that 
reason in Fig. 5 we show only upper limits on CIB flux 
corresponding to LCDM (black line) and PFI (grey line) 
models as representatives of these two group of models, 
respectively. 



is contributed by interactions with the 'dust' compo- 
nent of CIB at wavelengths A > 10/im (see Fig. 1). The 
optical depth increases with energy so rapidly that re- 
quiring the absorption-corrected spectrum to be concave 
is already sufficient to impose an interesting constraint 
on the CIB. Therefore the most robust upper limits on 
CIB at FIR are provided by the HEGRA data above 
10 TeV. Indeed, in this energy region the source spec- 
trum 'reconstruction factor' exp(r 77 ) should not exceed 
the exponential term of the observed spectrum of Mkn 501 
which in the energy region up to 24 TeV is well described 
as J obs (E) oc E~ 192 exp(-E/6.2TeV) (Aharonian et al. 
1999). The intrinsic (reconstructed) spectra of Mkn 501 
are shown in Fig. 3 assuming that the spectral shape of 
CIB is described by the LCDM model with 3 different 
scaling factors: SF = 1, 2.5, and 4. 

In order to obtain an upper limit on the CIB intensity 
around 50 /im, we require that the intrinsic flux at 17 TeV 
(the highest energy point in the HEGRA data set with 
adequate statistical significance; Aharonian et al. 1999) 
should not exceed the flux at 4 TeV. In the case of the 
SSC model this implies that for any reasonable combina- 
tion of model parameters, the Compton peak should not 
appear at energies beyond several TeV (in particular, due 
to the Klein-Nishina effect) . For the hadronic ( "7r°-decay" ) 
models (e.g. Dar & Laor 1997) this implies that the spec- 
trum of accelerated protons should not be flatter than the 
"nominal" E~ 2 spectrum. 

It is seen from Fig. 3 that the LCDM type spectrum 
of CIB with scaling factor up to SF — 2.5 do not con- 



Model 


S Fmax 


Upper Limit 


(nW m" 2 sr" 1 ) 






at 20 fj,m 


at 80 /mi 


LCDM 


3.4 


3.6 


21.1 


CHDM 


2.9 


2.9 


18.0 


PFC 


2.6 


5.9 


11.4 


PFI 


1.7 


6.3 


12.7 


ED 


1.15 


2.9 


14.2 



Table 1. The 68 % C.L upper limits on the scaling factor SF 
derived at FIR from 20 to 80 fim for different CIB models. 



Finally we note that although the extrapolation of the 
FIR upper limits shown in Fig. 5 do not contradict the 
COBE measurements at A > 100 /jm, the maximum en- 
ergy of about 17 TeV reported from Mkn 501 with high 
statistical significance, is not sufficient to provide model- 
independent information about the CIB at wavelengths 
beyond 80 [im. 

3.2. Limit on the CIB at NIR 

At shorter wavelengths, A < 10 /im, the constraints on the 
CIB come mainly from data below 10 TeV. For a power- 
law spectrum of the CIB, n(e) oc e - ' 3 , the optical depth is 
proportional to r 77 oc E" -1 (see e.g. Gould and Schreder 
1967). In the 2-10 /im side of the 'valley' the CIB, con- 
tributed mainly by the starlight, has a typical spectrum 
\F\ oc X s with S ~ — 1 (see e.g. Dwek et al. 1998; Primack 
et al. 1999), or ticib(£) oc e _1 . Therefore the shape of the 
7-ray spectrum at energies between 2 TeV and 10 TeV 
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remains essentially unchanged (r 77 (S) ~ constant), al- 
though the absolute absorption effect could be very large. 
This effect is clearly seen in Fig. 2. It makes rather dif- 
ficult to extract information about CIB based merely on 
the 7-ray data in the intermediate energy range of 7-rays 
between 2 and 10 TeV. 

Essentially more information about the CIB is con- 
tained in low energy, E < 1 — 2 TeV 7-rays, since in this 
energy region the absorption effect again becomes energy- 
dependent (see Fig. 2). For the redshift of Mkn 501, z = 
0.034, the optical depth t 77 at 500 GeV remains less than 
1 for any reasonable assumption about the CIB at NIR. At 
1-2 TeV the optical depth could be close or even exceed 1, 
if SF > 2. Therefore, a variation of the density of the CIB 
by the scaling factor SF from 1 to 4 could lead to dramatic 
(up to factor of 10) increase of the 'reconstructed' flux at 
1-2 TeV, while at energies below 500 GeV the impact is 
still not very large (less than factor of 2) . This implies that 
the slope of the 7-rays energy spectrum from Mkn 501 be- 
low 1 TeV is very sensitive to the level of the CIB flux at 
NIR. Therefore it can be used not only to constrain the 
CIB, but also to "measure" the CIB flux at wavelengths 
from 0.25 to few /xm, provided that a reliable intrinsic 
7-ray spectrum could be derived from multi-wavelength 
studies of the source. Note that the result is insensitive to 
the FIR density. 

The April 16, 1997 flare of Mkn 501 is particularly well 
suited for this task because (i) the appropriate distance to 
the source of about 170 Mpc - sufficiently large to provide 
non- negligible (e.g. measurable absorption at E < 1 TeV), 
and at the same time not too large for strong suppression 
of TeV emission; (ii) the availability of high quality X- 
ray and sub- TeV 7-ray data obtained simultaneously by 
BeppoSAX and CAT, and (iii) the general belief that 
the TeV radiation of this source has a synchrotron-self- 
Compton origin (SSC). In the SSC scenario the TeV 7- 
rays and synchrotron X-rays in a broad band from 0.1 keV 
to > 100 keV are produced by the same population of elec- 
trons accelerated in the relativistic jet with Dopplcr factor 
> 10 (for a review see Sambruna 1999). If the SSC model 
works for this source, it is possible to robustly predict the 
intrinsic 7-ray spectrum based on the multi-wavelength 
observations during a strong flare of Mkn 501. 

In this paper we do not intend to model the 
April 16, 1997 flare, but rather we use a few well known 
results from previous applications of the SSC model. 
Within this, for a realistic set of parameters character- 
izing the synchrotron X-ray and inverse-Compton (IC) 7- 
ray emitting jet, the 7-ray photons with energy < 1 TeV 
or so are produced in the Thomson regime, and therefore 
are described by the same spectral index as the X-rays do 
below the so-called break (see e.g. Tavecchio et al. 1998). 
The spectrum of X-rays of April 16, 1997 flare measured 
by BeppoSAX (Pian et al. 1998) had a photon index 1.5 
at 1 keV, and remained unusually hard up to very high 
energies with photon index 1.7 at 100 keV. Although for 



more definite conclusions one needs comprehensive mod- 
eling of the flare in a broad frequency band, we suppose 
that the sub- TeV 7-ray spectrum should repeat the spec- 
tral shape of relatively low energy X-rays, so we expect a 
photon index of intrinsic spectrum of sub- TeV 7-rays close 
to 1.5 (see e.g. Krawczynski et al. 1999 for more gen- 
eral discussion). Another characteristic feature of the 7- 
ray spectrum is its so-called 'Compton peak' which un- 
avoidably appears in any SSC or, more generally, external 
radiation Compton models (see e.g. Sikora et al. 1997). 

The spectrum of the IC 7-rays in a broad energy band, 
including the transition region from the Thomson regime 
to the Klein-Nishina regime (which is around 1 TeV here), 
has a rather complicated form. For convenience, following 
the recommendation of Sikora et al. (1997), we fit the 
7-ray spectrum, corrected for the intergalactic absorption 
with a given scaling factor SF, on the (log(vF(v)), log(E)) 
plane by a parabola with a slope at 600 GeV a^oo (photon 
index) . 

The relation between SF and a^oo computed by % 2 
minimization is shown in Fig. 4. These results imply that 
the photon index expected around 1.5 within SSC model 
for the 'Thomson' branch of the IC spectrum could be 
achieved for a scaling factor SF ranging between 0.5 and 
2.8 at 99% CL (LCDM model). This implies a CIB flux of 
ps5-35 nW m -2 sr -1 at 1 /an. The analogous calculations 
for each of the other four CIB models lead to the limits 
summarized in Tab. 2. Remarkably, the difference between 
the results derived for different models is always less than 
20-30 %. Thus the conclusion about the lower and upper 
limits of variation of the CIB intensity at 1 /im, allowed 
by the SSC model, is very robust. 

The above estimate of the CIB flux is based on a 
model-assumption that the observed sub- TeV 7-rays are 
produced within the framework of the SSC model as 
a result of inverse Compton scattering of highest en- 
ergy electrons responsible also for the synchrotron radi- 
ation of the jet. However, the strong relation between 
the scaling factor SF and the photon index of the in- 
trinsic sub- TeV 7-ray spectrum allows very robust con- 
straint on the CIB without a specific model assump- 
tion. Indeed from Fig. 4 follows that at large scaling fac- 
tor SF, the 'reconstructed' 7-ray spectrum becomes ex- 
tremely hard, e.g. with a6oo < 1 for SF w 5. Such 
spectrum seems to be unacceptably hard for any realistic 
scenario of 7-ray production connected with either pro- 
tons or electrons. So we may draw a conclusion that at 
99% C.L. the value ~ 60 nW/m 2 sr should be considered 
as an absolute upper limit on the CIB flux at 1 ^m. An as- 
sumption of a smaller distance to Mkn 501, e.g. adopting 
for the Hubble constant H =75 km.Mpc^.s -1 (instead of 
60 km.Mpc _1 .s _1 ), could soften this upper limit only by 
25 %. 

In Fig. 5 we compare the limits on the CIB intensity 
derived above with the direct measurements. Results ob- 
tained in this paper are shown by large horizontal bars 
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corresponding to the energy range of CIB photons which 
can interact with 600 GeV photons. The 99% CL limits 
correspond to the 5 different CIB models. The upper limit 
at 60 nW m~ 2 sr _1 assumes a differential source spectrum 
softer than E~ x , The allowed range of variation of the CIB 
between 5 and 35 nW m~ 2 sr _1 is derived assuming in- 
verse Compton origin of TeV radiation of Mkn 501. These 
results are in good agreement with recently reported fluxes 
of CIB at 2.2 /im and 3.5 (im based on the DIRBE mea- 
surements (Dwek & Arendt 1998, Gorjian et al. 1999), as 
well with the measurement of the background radiation 
from absolute photometry (Bernstein et al. 1999). 



Model 


Upper Limit 


SSC limits 




nW m~ 2 sr" 1 


nW m -2 sr -1 


LCDM 


48 


5 to 28 


CHDM 


60 


7 to 34 


ED 


51 


7 to 30 


PFC 


44 


5 to 26 


PFI 


57 


6 to 34 



Table 2. The upper limit on the CIB based on the assump- 
tion that the differential source spectrum of 7-rays at sub- TeV 
energies is softer than -E -1 , and the interval of CIB densities 
allowed by the SSC model of TeV radiation of Mkn 501. The 
99 % C.L. limits are relevant to the NIR wavelengths from 0.4 
to 3 fim, corresponding to the energy range of interaction for 
600 GeV photons. 



4. Impact of the CIB flux on the jet parameter- 
space of Mkn 501 

The analysis of the spectral shape and variability of the 
synchrotron and IC components of non-thermal radia- 
tion of TeV blazars within the framework of a single- 
zone SSC model may yield important constraints in the 
parameter-space of the X-ray and 7-ray production region 
(see e.g. TMG, Bednarek & Protheroe 1999). Among the 
key observables used in derivation of constraints on the 
jet parameters are the frequency vie and the apparent 
luminosity 2 vqL{vjc) of the Compton peak. Since the 
intergalactic absorption of 7-rays could significantly de- 
form the original information about both parameters (see 
Fig. 6), we study the impact of the CIB density in the 
parameter-space of the jet in Mkn 501. For calculations, 
we use the convenient analytical approach recently devel- 
oped by (TMG) for the homogeneous SSC model of TeV 
blazars. Detailed discussion of the SSC model is not the 

2 Below for convenience we use the apparent luminosity of the 
source which is defined as the luminosity of an isotropically 
emitting source at a given distance. The intrinsic luminosity 
of the blob relativistically moving towards the observer with 
Doppler factor 8 3> 1 is much (oc 5 4 ) smaller. 




0.5 1 1.5 2 2.5 



600 GeV differential photon index 

Fig. 4. Relation between the scaling factor SF and the photon 
index of the 7-ray spectrum after correction for absorption in 
the CIB, a6oo, obtained from the \ 2 fit for the LCDM model. 
The contour levels are at 10, 90 and 99 % C.L.; the scaling 
factor related to the upper limit is the maximum scaling factor 
compatible at 99 % C.L. with aaoo = 1 while the scaling fac- 
tors allowed by the SSC model lie between the minimum and 
the maximum scaling factors compatible at 99 % C.L. with 
CK600 = 1-5. 

aim of this paper but rather we show the impact of the 
CIB absorption on the derivation of these parameters. 

4-.1. Model parameters and observables 

The homogeneous SSC model describes the emission of 
electrons in a single blob with three basic model parame- 
ters: the radius R, magnetic field B, and Doppler factor of 
the bulk motion 6 of the blob. An important information 
about the ratio R/S is contained in the observed source 
variability time-scale t var , namely the radius of the source 
should not exceed R ma x — t var 'xcxS. Thus the constraints 
in the parameter-space for the given maximum size (de- 
fined through the variability time-scale) can be expressed 
in terms of two parameters on the (logB, log S) plane. 

A distinct feature of the SSC model is the character- 
istic SED of radiation with two pronounced, synchrotron 
and Compton bumps which, in the case of X-ray selected 
BL Lac population, appear in the X-ray and the TeV 7- 
ray bands respectively (see e.g. Ulrich et al. 1997). The 
observations of X-rays require a population of relativis- 
tic electrons with power-law spectrum broken at energy 
76. The specific values of the indices are determined by 
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Fig. 5. CIB spectral energy distribution. The results of this 
paper are represented by large bars (0.4-3 /xm, 99% CL) and 
are valid for all reasonable models or the SSC model (grey 
bars). Two limits are plotted in the range 20-80 /xm (68% CL) 
as the results depend on the CIB slope. See text for details. 
The DIRBE results are shown by black squares (Hauser et al. 
1998). Three grey squares correspond to CIB density extracted 
from WHAM H a survey and Leiden/Dwingeloo HI data at 
100, 140 and 240 /xm (Lagache et al. 1999). The grey square 
at 15 /xm comes from galaxies counting (Elbaz et al. 1998). 
The two grey dots show the fluxes of 2.2 and 3.5 /xm (Dwek 
& Arendt 1998, Gorjian et al. 1999). The lower limits corre- 
spond to the Hubble Deep Field galaxy counts (Pozzetti et al. 
1998), and the grey stars come from HDF results combined 
with ground-based spectrometry (Bernstein et al. 1999). The 
range of the CIB flux detected by FIRAS at very long wave- 
lengths is shown by dotted line (Fixsen et al. 1998). 

the spectral shape of the synchrotron radiation below and 
above the synchrotron peak. 

In the quiescent state, the synchrotron radiation of 
Mkn 501 is characterized by a synchrotron peak at v s ~ 
10 16 Hz, and spectral indices ol\ — 0.5 and a 2 = 1.75 
(TMG). During the April 16, 1997 flare the X-ray spec- 
trum was exceptionally hard with photon index less than 
2 at least up to 100 keV, indicating a dramatic shift of 
the synchrotron peak by at least two orders of magnitude 
(Pian et al. 1998). Due to the lack of statistics, both the 
exact position for the synchrotron peak position and the 
synchrotron luminosity during the April 16, flare are not 
well defined. 

The synchrotron spectrum of Mkn 501 in the high 
state is represented in the form of broken power-law with 
spectral indices ot\ = 0.5 and a.% = 1 below and above 
the energy hv b = 21.5 keV {v h = 5.2 ± 0.3 x 10 18 Hz (see 
Fig. 6). Additionally, we assume a high energy cutoff at 
300 keV, which naturally could be attributed to the cut- 
off in the acceleration spectrum of electrons. The index 
c*2 as well as the position of the cutoff are rather qualita- 
tive, but fortunately the final conclusions do not depend 




v (Hz) io" io 1a 10" io 20 v (Hz) io 35 io 2E io !7 io ! 



Fig. 6. Spectral energy distribution of the April 16, 1997 flare 
of Mkn 501. The X-ray data are from BeppoSAX observations 
(Pian et al. 1998). The 7-ray fluxes measured by CAT are 
shown by filled squares. The HEGRA 'time-averaged' spectral 
points obtained during the entire 1997 outburst of Mkn 501 
(Aharonian et al. 1999) are also shown (open squares). Statis- 
tical errors only are drawn. The filled dots correspond to the 
CIB absorption-corrected fluxes computed for LCDM model 
(SF — 2.5). The 3 a are used to model the spectra. 

strongly of their exact values. The total luminosity is then 
v s L(v s )=2.b±0.l 10 45 erg.s" 1 . 

According to the TMG description, the IC peak should 
be represented by a broken power-law with the same spec- 
tral index ot\ = 0.5 before the maximum and the spectral 
index a 2 — 1.5 after, as shown in Fig. 6. The position of 
the maximum, and therefore the total energy depend on 
the assumption about the CIB absorption. Tab. 3 presents 
these observable parameters computed with no absorption 
and with the LCDM and PFI models, using SF=1 or 2.5. 

It is interesting to note that the correction of the ob- 
served 7-ray spectrum for the intergalactic absorption, as- 
suming typically a CIB flux at the level of 25 nW ~ 2 sr _1 
at A ~ 1 /xm, results in a symmetric SED with a distinct 
maximum at 2 TeV, and power-law type behavior with 
spectral indices 0.5 and 1.5 below and above the peak, re- 
spectively, although we constrained only ct\ — 0.5. This is 
a feature predicted by the SSC model for the IC radiation 
provided that the synchrotron X-radiation is described by 
a broken power-law with spectral indices a\ = 0.5 and 
a 2 = 1 (TMG, Krawczynski et al. 1999). 

4-2. The constraint regions on the (log_B, logd) plane 

TMG have suggested three independent constraints on al- 
lowed regions in the (log B, log <5) plane based on (i) the 
positions of the synchrotron and Compton peak frequen- 
cies (area A), (ii) the synchrotron and Compton peak 
luminosities (area B), and (iii) the equilibrium between 
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Model v IC i>icF{vw) E IC , na x vicL{v IC ) 

no absorpt. 2.1±0.2 3.9±0.1 0.9T0.1 1.3±0.1 

LCDM xl 3.6±0.2 6.0±0.2 1.5±0.1 2.1±0.1 

LCDM x2.5 5.3±0.3 13.1±0.3 2.2±0.1 4.5±0.1 

PFI xl 4.5±0.2 9.4±0.2 1.9T0.1 3.2±0.1 

PFI x2.5 7.7±0.5 40.0±0.1 3.2±0.2 13.9±0.4 



total energies of the synchrotron and Compton peaks as 
described above. Also, it is interesting to note that the 
ignorance of the CIB absorption of 7-rays would lead to 
a conclusion that the radiative cooling of electrons in the 
jet is well dominated by synchrotron losses. This, however, 
could not be true, since after the correction of 7-ray fluxes 
for the intergalactic absorption, the IC luminosity in fact 
could be as high as (or even exceed) the synchrotron lu- 
minosity (see Fig. 6). 

A relation between B and <5 arises if one assumes 
that the break in the electron spectrum at 7^ is deter- 
mined from the equilibrium between the cooling and es- 
cape from the source (obviously, there could be other rea- 
sons for the break in the electron spectrum, e.g. connected 
with the character of the acceleration mechanism). This 
model assumption leads to the following constraints in 
the Compton-cooling dominated regime (eq. 32 and 34 in 
TMG) 4 : 

(CI) B8 9 10 > 37.(^ s ) 45 ) 2 Cfc &~>7c,26 (3) 

and synchrotron-cooling dominated regime: 

(<72) B > 0.18 /£/ c 2 (W,h ^c,2 6 )" 1/2 (4) 

Here /3 esc is the electron escape velocity in units of 
the speed of light, i.e. a parameter which describes the 
energy-independent escape time as t esc = [3 esc R/c. Fol- 
lowing TMG, we allow a change of this rather uncertain 
parameter within limits from 1/3 to 1. 

To avoid strong absorption of TeV photons inside the 
blob due to pair-production of 7-rays interacting with op- 
tical photons, a minimum value for 5 can be computed ; 
it is used to check the validity of the area defined by the 
intersection between (A), (B) and (C) regions. For the 
Compton peak luminosity, we obtain 

(£>) 5 > 6.6 (^l)^ V v Ica6 ^ 19 ) ° (5) 

The impact of the intergalactic absorption of 7-rays on 
the constraints described by the regions (B), (C) and (D) 
is less dramatic than on the region (A). Nevertheless, as 
can be seen from Fig. 7, even for these regions the ef- 
fect is not negligible, and should be taken into account 
in any realistic attempt to constrain self-consistcntly the 
parameter-space in future studies based on detailed mod- 
cling of temporal and spectral characteristics of radiation. 
The scaling factor SF which determines effectively the 
flux of CIB at NIR, and therefore defines the level of dis- 
tortion of 7-radiation in the sensitive region of the Comp- 
ton peak, should be considered as a new free parameter in 
such studies. 

Although the detailed modeling of radiation of 
Mkn 501 is beyond the framework of this paper, our study 

4 There is a mis-print in the TMG paper. The exponent of 5 
in Eq. 34 should be (6/(1 - ati) - 3) instead of (6/(1 - ai)). 



Table 3. The parameter Vic (in units 10 26 Hz), vwF{vic) 
(in units 10 -9 erg cm~ 2 s _1 ), Ejcmax (in units 10 12 eV) and 
vi C L(yic) (in units 10 45 erg s _1 ) for two CIB models (LCDM 
and PFI) and two values for the scaling factor SF (1 and 2.5). 
The first line of the table corresponds to the 7-ray spectrum 
without intergalactic absorption. 



radiative cooling and escape of electrons (area C). It is 
important to note that the convenient approximate an- 
alytical formulae derived by TMG are applicable to the 
treatment of the IC scattering relativistic (Klein-Nishina) 
regime required here. 

The first condition leads to the following simple rela- 
tion between B and S (eq. 16 in TMG) 3 : 



(A) B s£Ki. 5 Jp»- (1) 



where B is in Gauss, Sio = 5/10, v s ,i9 = ^ s /10 19 Hz, and 

view = vic/10 26 ~Rz- 

The IC peak is very sensitive to the CIB absorption. 
Typically, the peak position shift by a factor 3 (up to 10) 
towards higher frequencies for the absorption-corrected 
spectrum. The total energy follow a similar evolution. 
This implies that the ignorance of the intergalactic 7- 
ray absorption would overestimate the B/S ratio by a fac- 
tor of up to 10 or even more. This is seen in Fig. 7 from 
comparison of the regions (A) on the top and bottom pan- 
els. 

The comparison of the observed luminosities in the 
synchrotron and the IC peaks gives the second relation 
between the magnetic field and the Doppler factor (eq. 22 
in TMG) 

(B) B5\i > 0.5 K, 19 „ JOia6 )-V i ^ s) \ i " C,,(2) 

where {v s L(v s ))±5 = {v s L(v s )) /I0 ib erg s" 1 , 
(vicL{v IC ))ia = (vicL(v ic ))/10^ erg s^ 1 are the syn- 
chrotron and IC peak apparent luminosities and t var ^ = 
t var /l h is the source variability time-scale. Here t vart h=10 
is assumed. 

The region (B) shown in Fig. 7 is based on Eq. (2) with 
allowed ranges of uncertainties (3 a) in the positions and 

3 Since we do not have here 7mai(=300 keV) S> 7b(=100 
keV), the expression of g(ai, 012) must be modified. g(ai, Q2) = 
exp(l/(ai - 1) + 1/(2 * (a 2 - ai)) * (1 - v / j b /'yma X ) 
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shows that the inclusion of the intergalactic absorption in 
the treatment results in a more consistent picture. In par- 
ticular it moves the region (A) into the regions (B) and 
(C), and thus overcomes the incompatibility of different 
constraints on the principal model parameters. It is seen 
from Fig. 7 that after the correction of 7-ray fluxes for 
the intergalactic absorption, there is a trend for appear- 
ance of a compact common area for all three (A), (B), and 
(C) regions in the (log B,log S) plane with magnetic field 
B ~ 0.05 G, and Doppler factor 5 ~ 15. 
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Fig. 7. Parameter-space for the April 16, 1997 flare of 
Mkn 501 (with t var ~Wh) and SF— (i.e. no intergalactic 
absorption) or SF = 2.5 for the LCDM model. See text for 
details on areas (A), (B), (C) and (D). 



If the SSC model is confirmed in future studies, this 
would be also a strong evidence for the high CIB at 
NIR as is shown in Fig. 5, and used for the reconstruc- 
tion of the source spectrum of the April 16, 1997 flare 
in Fig. 6. Such a high flux of CIB makes significantly 
steeper (Acv ~ 0.5) the spectrum of sub-TeV radiation 
even from a relatively nearby source like Mkn 501. Ob- 
viously, the impact to the 7-ray emission from farther 
sources could be much stronger. In Fig. 8, the intergalac- 
tic absorption factor, exp(— t 77 ) is shown for Mkn 501 and 
for the BL Lac object PKS 2155-304 (z=0.116), which re- 
cently has been reported by the Durham group as a 7- 
ray emitter above 300 GeV (Chadwick et al. 1999). It 
is seen that for a high CIB flux with the scaling fac- 
tor SF — 2.5 the 7-ray spectrum should suffer dramatic 
steepening (Aa ~ 2), before reaching the observer. Thus, 
even for a quite flat source spectrum, e.g. with photon 
index 1.5, we should expect very steep 7-ray spectrum 
from this source, (dN/dE cx E~ 3 5 ). Interestingly, at en- 
ergies above 2 TeV the absorption becomes almost energy 
independent, therefore the observed spectrum essentially 
repeats the spectral shape of the source spectrum. How- 
ever, the detection of > 2 TeV radiation of 7-rays from 
PKS 2155-304 would be very difficult because of strong, 
by a factor of 100, suppression of 7-ray fluxes in that en- 
ergy regime. 




1 Energy (TeV) 10 

Fig. 8. Intergalactic absorption factor for Mkn 501 (z=0.034) 
and PKS 2155-304 (z=0.116) calculated for the LCDM model 
assuming two different scaling factors, SF — 1, and SF = 2.5. 



5. Summary 

The accurate spectromctric observations of the 7-ray spec- 
trum of Mkn 501 during the remarkable April 16, 1997 
flare by the low-threshold imaging atmospheric Cherenkov 
telescope CAT in the energy region from 300 GeV to 
10 TeV is used for extraction of important information 
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about the CIB at wavelengths > 0.4 fim. The interpre- 
tation of the spectrum of sub-TeV 7-rays, together with 
simultaneously obtained X-ray data of BeppoSAX, re- 
quires, within the one-zone SSC model, rather high NIR 
background at a level close to 20nWm~ 2 sr~ 1 at 1 /zm. 
Such high flux of CIB implies an essential distortion of 
the shape of the initial (source) spectrum of 7-rays from 
Mkn 501 not only at multi-TeV, but also at sub-TeV ener- 
gies. The 'reconstructed' intrinsic 7-ray spectrum shows a 
distinct peak in the Spectral Energy Distribution around 
2 TeV with a flux by a factor of 3 higher than the mea- 
sured flux. Moreover, the energy spectrum of gamma radi- 
ation from both sides of the peak has power-law behavior 
with spectral indices a ~ 0.5 below 2 TeV, and a ~ 1.5 
above 2 TeV, which perfectly agrees with predictions of 
the SSC model. We have shown that the intergalactic ab- 
sorption has non-negligible impact on the construction of 
self-consistent SSC parameters. 

And finally, we argue that the CAT 7-ray data alone 
allow rather robust upper limits on the CIB, XF\ < 
60 nWm^sr 1 at 1 /zm, taking into account that for 
any reasonable scenario of 7-ray production the differen- 
tial intrinsic spectrum of 7-rays hardly could be flatter 
than diV/dSoc E~ x . 
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